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(57) Abstract: The invention relates to measuring instruments, preferably of the kind measuring absorbances, in an object, of elec- 
tromagnetic radiation in at least two spectral ranges, such as IR instruments, and DXR, meaning Dual X-ray instruments, and more 
specifically to the determination of properties of food or feed, such as the fat content in milk or meat The invention relates in par- 
ticular to a method of providing a correction for a slave instrument of the kind measuring properties of an object by exposing the 
object to electromagnetic radiation, in particular X-rays, in at least two spectral ranges and obtaining one or more object responses 
thereto. The responses obtained being preferably based on detecting attenuation and/or reflection and/or scatter of the electromag- 
netic radiation in/from the object by use of one or more detectors and are obtained in a form where they express properties of the 
object either directly or via a transformation. 
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Title 

Method and means for correcting measuring instruments. 
5 Technical field 

The present invention relates to measuring instruments, preferably of the kind measuring 
absorbances, in an object, of electromagnetic radiation in at least two spectral ranges, 
such as IR instruments, and DXR, meaning Dual X-ray instruments, and more 
10 specifically to the determination of properties of food or feed, such as the fat content in 
milk or meat. 

Background art 

Spectral instruments measuring e.g. infrared absorbances at several wavelengths in 
15 order to determine contents of specific components in a liquid such as milk are well 
known. Also X-ray analysis for determining the fat content of meat has been known for 
several years. 

Typically such instruments apply regression analysis and multivariate calibration. Such 
20 analysis is known from e.g. the applicant's own PCT application No. WO 95/16201 
disclosing the determination of extraneous water in milk samples using regression 
analysis and multivariate calibration. Further, the applicant's PCT application No. WO 
98/43070 discloses measurement of acetone in milk using IR spectroscopy and 
multivariate calibration. The transfer of calibrations from one instrument to another has 
25 been discussed in US patent No. US 5,459,677 disclosing a "Calibration transfer for 

analytical instruments" and US patent No. US 5,559,728 disclosing "Calibration transfer 
for second order analytical instruments" and in the applicants US patent No. US 
5,933,792 "Method of standardizing a spectrometer* 

30 The applicant's WO 93/06460 discloses an infrared attenuation measuring system, 
including data processing based on multivariate calibration techniques, and the 
applicants US 5,252,829 discloses a determination of urea in milk with improved 
accuracy using at least part of an infrared spectrum. 

35 As disclosed in WO 01/29557 the properties a medium of food or feed, such as the fat 
content of meat, may be determined by use of dual X-ray absorptiometry, the medium 
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being a raw material of food or feed, a product or intermediary product of food or feed, 
or a batch, sample or section of the same, the method comprising - scanning 
substantially all of the medium by X-ray beams having at least two energy levels, 
including a low level and a high level, - detecting the X-ray beams having passed 
5 through the medium for a plurality of areas (pixels) of the medium, - for each area 
calculating a value, A| 0W , representing the absorbance in the area of the medium at the 
low energy level, - for each area calculating a value, Aw gh representing the absorbance 
in the area of the medium at the high energy level, further comprising for each area 
generating a plurality of values being products of the type A {ow n * A high m wherein n and 

10 m are positive and/or negative integers or zero, and predicting the properties of the 

medium in this area by applying a multivariate calibration model to the plurality of values, 
wherein the calibration model defines relations between the plurality of values and 
properties of the medium. The advantage over the prior art is a more accurate 
determination of the properties, such as the fat content in the medium. The accuracy is 

15 specifically improved over the prior art when measuring layers of varying thickness. A 
further advantage is due to the fact that using the described method almost the whole 
product is measured instead of a sample thereof. Generally, extraction of a sample from- 
an inhomogeneous medium will introduce an error, because the sample may not be truly 
representative. 

20 

Preferably the calibration model is obtained by use of a multivariate regression method 
being included in the group comprising Principal Component Regression (PCR), 
Multiple Linear Regression (MLR), Partial Least Squares (PLS) regression, and 
Artificial Neural Networks (ANN). 

25 

A problem related to the prior art: 

It is well known, that when a number of measurement instruments measure the same 
sample, each instrument will generally produce an instrument specific signal if no specific 
actions are taken to ensure that the signals produced by the instruments are identical for 
30 an identical sample. It is equally well known that it is desirable to. be able to manufacture 
measuring instruments, which generate the same signal when exposed to the same 
sample. 

Calibration of an instrument may remedy the problem. However the multivariate 
35 calibration, which is applied for the DXR analysis - as described above and in the 

published WO 01/29557 - is a delicate matter requiring a number of known reference 
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samples, which typically have to be analyzed by an officially recognized reference 
method. The provision of such analysis results of the delicate perishable calibration 
samples consisting of various mixed samples of meat and fat which have to be handled 
very carefully preferably as frozen items - are time consuming, tedious and expensive. 
Further the calculations required for providing the calibration are time consuming too and 
therefore, expensive. These drawbacks are emphasized when a large number of 
instruments must be calibrated. Furthermore, when such calibrations must be performed 
often, such as regularly due to drifting in the instruments, the calibration method may 
strongly influence the usefulness of the instruments in a negative manner. 



Thus, an object of the present invention is to provide a method, which seeks to avoid the 
above-mentioned drawbacks. Accordingly it is an object to disclose a method and 
instruments enabling a single highly sophisticated calibration developed on a master 
instrument to be applied to all other instruments in a series of similar instruments. 



More specifically it is an object to provide a method of adjustment or correction for a 
series of instruments in such a manner that they can use the same calibration. 

The present invention is specifically useful to measurements on meat performed by use 
of Dual X-ray equipment designed for measuring fat and areal density in meat, as well as 
for detecting foreign bodies in a meat sample. Such a measurement should, in order to 
obtain an acceptable accuracy, detect X-ray attenuation at at least two X-ray energies. 
According to the particular aspect of the present invention an X-ray equipment 
comprising two X-ray sources and two X-ray detectors is used for measuring the 
absorbances. Measurements performed with such equipment have shown to be 
extremely delicate as the amount of X-rays absorbed by adipose (fat) and muscle tissue 
only differs slightly, thereby demanding extreme care in calibrating the instrument. 

Thus, another object of the present invention is to provide a method, which ease set-up 
of an Dual-Xrray instrument so that it is capable of producing accurate measurement of 
fat and areal density in meat. 

A further object is to allow for less restrictive instrument specifications, enabling the use 
of cheaper components having rather coarse tolerances, thereby reducing the total cost 
of the instrument. 
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Disclosure of the invention 

In a first aspect of the present invention and in accordance with the objects of the 
invention, a method of providing a correcting for a slave instrument is suggested. The 
5 slave instrument is preferably of the kind measuring properties of an object by exposing 
the object to electromagnetic radiation, in particular X-rays, in at least two spectral 
ranges and obtaining one or more object responses thereto. The responses obtained 
being preferably based on detecting attenuation and/or reflection and/or scatter of the 
electromagnetic radiation in/from the object by use of one or more detectors and are 
10 obtained in a form where they express properties of the object either directly or via a 
transformation. The suggested method of correcting comprises: 

obtaining, for a plurality of stable objects, a set of responses comprising one or 
more pair of related responses (Qiow s and Qwgh 5 ) representing measurements in 
the at least two spectral ranges performed with the slave instrument and a set of 
15 responses, comprising one or more pair of related responses (Qiow m and Qhigh™) 

representing measurements in the at least two spectral ranges performed with a 
master instrument, 

to each pair of related responses (Qiow s and Q h j g h s ) of the slave instrument 
corresponds a pair of related responses (Qiow m and QhigO of the master 
20 instrument, 

and to each element in each pair of responses (CW and Qwgh 5 ) of the 
slave instrument corresponds an element in the corresponding pair of 
responses (Qiow 01 and QngtF) of the master instrument; 

- determining based on the sets of responses a correcting function being a • 
25 functional relationship between a ratio of related responses of the master 

instrument and a sum of a plurality of terms, each term being a product of a 
correcting coefficient (Bj) and powers of related responses (Q tow s and Qhigh 8 ) of 
the slave instrument, wherein each response being raised to a power being a 
positive or negative real number, or zero, thereby determining a first set of 

30 v correcting coefficients (B 0 ; B^- B 2 ) being multiplied on each of the terms; 

and 

- storing the first set of correcting coefficients (B 0 ; Bi; B 2 ) in memory means 

included in or adapted for communication with data processing means included in 
or adapted for communication with the slave instrument 



35 
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Thanks to the provision of a correcting function, each instrument may be brought in a 
condition where each of them produces similar, such as identical, responses on identical 
objects. Thus, by utilizing the correcting functions the responses produced by slave 
instruments can be transformed into so-called standardized responses (or corrected 
5 responses) which when subjected to a calibration function depending on such 
standardized responses may provide the desired knowledge about the physical 
properties of a measured object. 

Thus, while calibrations typically are tedious and expensive due to required chemical 
analysis of each object used in calibration and due to the high number of calibration 
objects needed for a satisfying calibration, the correcting function may be obtained by 
use of a set of standardizing objects and a corresponding set of master responses 
preferably stored in any suitable medium, preferably a data memory, such as a disc 
delivered together with the set of standardizing objects or alternatively, a report written 
on paper. 

Thereby, the most common job to be performed in connection with the present invention- 
is to achieve the same responses from the slave instrument, as would have been 
provided by the master instrument measuring the same object. This advantage is 
20 believed to be achieved by the present invention. 

Thus, a calibration model will thereby be transferable between all standardized 
instruments, and the measurements, such as the X-ray determination of fat in meat, will 
be much easier to handle. The present method for correcting a slave instrument, such as 
25 a Dual-X-ray equipment, using a limited number of stable standardizing objects, is 

therefore considered very useful, and therefore also the method of obtaining a correcting 
function used in performing the correction. 

Thanks to the saved set of correcting coefficients (B 0 ; B 2 ) the slave instrument 

30 inclusive the data processing means and stored programs will able to correct successive 
measurements on unknown objects providing substantially correct measurements results 
based on a calibration elaborated for a master instrument. This will be explained further 
in the detailed description. 



35 



The term response as used in this application means is always related to a detected 
signal and generally it will be related to at least two detected signals. Accordingly the 
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response may be considered to be a mathematical transformation of a number of signals 
resulting from a detector. Such signals may typically be digital signals provided by an 
analogue to digital converter converting electrical signals, provided by the detectors, into 
digital representations of the signals. Thus, a response, for instance the intensity (I), is in 

5 this connection typically generated from a detector signal, typically being a voltage, 
current or digital representation, by applying a mathematical relation to the detector 
signal, such as l=f(U) where f is a mathematical function and U is the voltage provided by 
a detector. In that sense all the responses considered are typically resulting from 
transformations of signals into responses. However, it is contemplated that the invention 

10 is applicable also in embodiments where the signals from the detectors, for instance a 
digital representation (e.g. as binary numbers), are used as responses as these signals, 
of course, express properties of the object . 

Depending on the mode of employment of the invention, such responses are either used 
15 as they are provided, for instance by mathematical relations, or are transformed into 
transformed responses. Typical examples of responses are intensity, which may be 
considered as a response used as it is provided, transmittance through an object being - 
derived from measuring the intensity with and without the object, and therefore preferably 
being considered as a derived response, absorbance being derived from the 
20 transmittance, reflectance, which also may be considered as a response used as it is 
provided, and a Kubelka-Munk transformation being applied to the reflectance. 

It should be noted that the designations "high" and "low" are used in general for 
designating two values where one of these is higher relatively to the other value. 
25 Furthermore, two or more separate sources may be utilized in connection with the 

present invention as well a combination of one source and two filters emitting a low and a 
high energy beam. 

The method according to the present invention comprises preferably the step of initially 
30 at a manufactures site measuring the plurality of stable objects on a 

master instrument, thereby obtaining the set of responses representing 
measurements performed with the master instrument (Qiow m and QhigO, 
storing the set of responses (CW m and Q hIa „ m ) as a set of constant values 
in memory means, which is accessible from a slave instrument, when 
35 measuring the corresponding stable objects on a slave instrument in order 

to carry out a method of correcting according to the present invention. 
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Thereby, the responses needed for providing the correcting function is easy accessible 
and the workload connected with the method may be reduced. 

Typically and preferably, the set of responses measured by the master instrument is 
5 stored in memory means included in or adapted for communication with data processing 
means included in or adapted for communication with the slave instrument 

The so-called "standardization" or correction is preferably based on a set of stable items 
which initially have been measured on a master instrument and all responses have been 

10 recorded on a recording medium such as a disc or ROM, which in the future follows the 
specific set of stable objects. In view of the fact that such a set of stable objects typically 
is expensive a specific set of stable object could be used for many slave instruments, 
e.g. applied once a year during a maintenance visit by a service technician or when 
major changes are made to the instrument (such as a change of radiation source or a 

15 detector). In preferred embodiments such standardized responses and/or any specific 
calibration may be accessible through computer means of well-known art, e.g. through 
the Internet on a pay per use principle. 

Typically and preferably, the determination of the correcting function being based on a 
20 regression method, which has proven to be a very efficient manner to obtain the 

correcting function. In preferred embodiments, the regression method is selected from 
the group consisting of principal component regression, multiple linear regression, partial 
least squares regression, and artificial neural networks. The partial least square method 
.' has proven to be especially useful. 

25 . 

It is generally preferred, that the correcting function comprises a plurality of terms of the 
following form Q tow n1 * Qhigh" 1 wherein n1 and m2 are real numbers and/or integers, and 
n1 is positive. In accordance with preferred and practical very useful embodiments it is 
preferred that the correcting function comprises at least three of the following terms: CW 
30 , Q h igh, Qiow 2 , Qhigh 2 and Q low / Q Wgh . In order to achieve a higher accuracy in the 

correcting, more terms may preferably be added to the correcting function such that the 
correcting function comprises at least three of the following terms: 



35 



Qtow* Qhigh; Qtow 2 *Qhlgh; Qlow * Qhigh 2 ! Qlow 2 /Qhigh; Qlow / Qhigh 2 ; Q|ow 2 /Qhigh 2 ; Qlow 2 / Qhigh 2 . 
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In particular preferred embodiments of the invention, the correcting function is of the 
form: 

^ = B,Q s lovl +B 2 Q s high +B3QL 2 +B 4 Qy gb 2 + 6^^ +B 6 Q? 0W 2 QLg h +B 7 Q? 0W Q s high 2 

Qwgh 



Chigh 



high 



high 



QS 



low 



Qhigh _ 



+ B„ 



The above-defined method can be adequate for many cases. However, in a generalized 
version of the method the above-mentioned correcting function is accompanied by a 
further correcting function being determined based on the sets of responses and being a 
functional relationship between responses of the slave instrument (Qiow s or Q h i 9 h s ) and 
related responses (Q low m or Q h i 9 h m ) of the master instrument. Thus, the method 
preferably comprises determination of a second set of correcting coefficients (a; p). 

Experience has shown that in some cases the further correcting function may improve 
the correction. 

Preferably and typically, the further correcting function is a functional relationship 
between a high energy response of the slave instrument (Qwgh S ) and the related high 
energy response (Qhigh 1 ") of the master instrument. Furthermore, it is preferred that the 
further correcting function is determined by use of univariate linear regression. 

In accordance with preferred embodiments of the present invention the further correcting 
function is preferably of the form Q h igh m = ct-Qtugh 5 + fl- 
it is preferred that the set of responses for the master instrument and the set of 
responses for the slave instrument each comprises one pair of related responses for 
each stable object comprised in the plurality of stable objects. 

In many practical implementations of the method the related responses are 
advantageously obtained based on measuring on objects being conveyed. 

In cases where the detector or detectors used for providing the responses is/are 
sufficiently stable in time in the sense that it is not necessary to take measures to 
eliminate for instance detector drift in order to maintain the overall accuracy of the 
measurement, it might not be necessary to apply the method to any transformed 
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responses. In such cases it is preferred that each of the. responses (Q) is an intensity (I), 
if necessary corrected with respect to a variable parameter, e.g. dark current of the 
detectors. Thus, it might be preferred to correct the raw intensity with respect to the dark 
current in order to increase the stability of the responses by compensating for instability 
in the detector. 

Such intensities are especially useful in situations where the intensities vary linearly or 
substantially linearly with physical properties reflected, which in particular is the situation 
where the measured absorption characteristics of an object varies over a narrow interval. 

To compensate for further instabilities, e.g. due to radiation source instability, it is 
typically preferred that each of the responses is a transmittance (T) being derived from 
intensity as a ratio between intensity resulting from measuring on an object and 
reference intensity. Of course, such transmittances may also be preferred in general. 

In preferred embodiments a linearization of transmittance is applied. In such and of 
course other embodiments as well, it is preferred that each of the responses is an 
absorbance being defined as the negative logarithm to a transmittance (A=-log(T)) such 
as logarithm base 10, the natural logarithm, or any other logarithmic function. 

In particular preferred embodiments of the present invention, the responses for both the 
master and the slave instruments are absorbances being determined by calculating 



wherein the intensities (I) are obtained in a measuring region of the master instrument 

respectively the slave instrument by: 

exposing the object in the measuring region to low and high X-ray energies and 
detecting with detectors the intensities l S a mp ie(low) and lsan,pte(high) respectively 
detecting the intensities l da *(low) and l^high) from said detectors when no 
radiation reaches them; 

and 

exposing said detectors to the low and high X-ray energies signals when no 
object is present in the measuring region and detecting lair(low) and Whigh). 
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In other preferred embodiments, each of the responses is selected to be a reflectance 
(R) expressing the reflectance from the surface of the object. The measured reflectance 
being useful in situations where the object to measured has such a nature that the 
properties to be measured is expressed by the reflectance of the object. It is 
5 contemplated that the reflectance covers surface reflectance of the object as well as 
reflectance in general of the object. 

In preferred embodiments using reflectance it is preferred that the reflectance is 
linearized. In such embodiments the reflectance (R) is linearized, preferably by using the 
10 Kubelka-Munk transformation (K/S=(1-R)/2R). 

The present invention relates in a second aspect to a method of correcting responses 
representing measurements performed with a slave instrument, said method comprising 
for an object 

15 - determining based on measurements with the slave instrument a pair of related 
responses (Qiow s and Q h igh s ), 

determining the ratio [CW Qwgh] 00 " by a correcting function being a functional 
relationship between a ratio of related responses of the master instrument and a 
sum of a plurality of terms, each term being a product of a correcting coefficient 
20 (Bj) and powers of related responses (Q Iow s and Q Wg h s ) of the slave instrument 

wherein each response being raised to a power being a positive or negative real 
number, or zero, 

providing Q h fgh corr either by assuming that Q h igh co ' r is substantially equal to Q Wgh s or 
by use of a further correcting function correlating Qhigh conr with Q h i g h s ; 

25 and 

calculating 0*™ as Qhigh 00 ^ * [CW Q hW J C0,r \ 
thereby providing a set of corrected responses. 

Thus, in accordance with the invention, the method of correcting responses comprises 
30 preferably three operations, namely determination of a ratio of responses, determining 
one of the responses of the ratio of responses and finally multiplying said one of the 
responses on the ratio determined. It should be noted, that the assumption that Qhi g h corr 
is substantially equal to Q W gh s can be construed as being a further correcting function (of 
the form Qhigh 00 "* = Qwgh 8 ) also, but is not termed so in order to ease the following 
35 discussion only. 
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When a further correcting function is utilised, this further correcting function is preferably 
of the form Qhigh 00 * = a-Q high s + 0. It is contemplated that one or both of the coefficients (a, 
0) and especially 0 can be determined or selected to be equal to zero. 



10 



It is in general preferred that the correcting function comprises terms of the following 
form Q, ow n1 * Qhigh" 11 wherein n1 and m2 are real numbers and/or integers, and wherein n1 
is positive. In accordance with preferred and practical very useful embodiments of the 
present invention it is preferred that the correcting function comprises the following 
terms: Q low , Qhigh, Qiow 2 , Qhigh 2 and Qiow / Q W g h . In order to gain a higher accuracy in the 
correcting the correcting function may preferably comprise more terms such that the . 
correcting function preferably comprises the following terms: Q Iow . Q hig h; Qiow 2 *Qhigh; Qiow 

* Qhigh 2 ! Qlow 2 /QhIgh; Qiow / Qhigh 2 ; Qlow 2 /Qhlgh 2 ; Qiow 2 / Qhigh 2 - 



In particular preferred embodiments of the second aspect of the present invention the 
15 correcting function is of the form: 



f \corr 
Qiow. 

^ Qhigh j 



= B,Qf ow +B 2 Q^ +B 3 Qf ow 2 +B 4 Q^ 2 +B 5 Q I s ow Q^ gb +B 6 Q? 0W 2 C& gh +B 7 Q? ow Q 



S * 

high 



o 5 o s o s 



high 



Qhis 



Vhigh 



low 



high 



+ B ft 



Wherein the B's preferably are constants being in general real numbers. 

20 In preferred embodiments and as disclosed in relation to the first aspect of the present 
invention each of the responses (Q) is in some situations preferred to be an intensity (I), 
if necessary corrected with respect to the dark current of the detectors. Such intensity 
may advantageously be transformed into a transmittance (T) being derived from intensity 
as a ratio between intensity resulting from measuring on an object and a reference 

25 intensity. 



30 



In order to, for instance, linearize the responses, it might be preferred that each of 
responses is an absorbance being defined as the negative logarithm to a transmittance 
(A=-log(T)) such as logarithm base 10, the natural logarithm, or any other logarithmic 
function. 



In particular preferred embodiments the responses are absorbances being determined 
by calculating 
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risampleC^W)- 1 ^ 

tow ~ 8,0 L uc^-i^oow) 



and 

da* (low) J 
- I '"'sample (high) -Id^ (high) 

A "*~~ g, °L ^(high)-!^^) J 
wherein the intensities (I) are obtained in a measuring region of the slave instrument by: 
exposing an object in the measuring region to low and high X-ray energies and 
detecting with detectors the intensities IsampieOow) and Isampie(high) respectively 
detecting with the detectors the intensities Inflow) and l d a*(high) from said 
detectors when no radiation reaches them; 

and 

exposing said detectors to the low and high X-ray energies signals when no 
object is present in the measuring region and detecting l alr (low) and l a!r (high). 

Also in this aspect of the present invention it may be preferred that each of the 
responses is a reflectance (R) expressing the reflectance from the surface of the object 
and the reflectance (R) may preferably and advantageously be linearized, preferably by _ 
using the Kubelka-Munk transformation (K/S=(1-R)/2R). 

The correcting function and the further correcting function utilized in the second aspect of 
the present invention are preferably determined by the method according to the first 
aspect of the present invention. 

The present invention relates in a third aspect to a method of determining a physical 
quantity for an object by a slave instrument. In this aspect the method comprises 
preferably 

determining for the object corrected high and low energy responses (Oh** 0 * and 
Qtow 00 ") by utilizing the method according to the second aspect of the present 
invention, 

determining the physical quantity by applying on said corrected responses a 
calibrated functional relationship between (W™ and Q im T n and a physical 
quantity. 

In accordance with the third aspect, it is preferred that the calibrated functional 
relationship being a functional relationship between a physical quantity (such as the 
content of a specific constituent e.g. fat and meat), and a sum of a plurality of terms, 
each term being a product of a calibration coefficient (Bj) and powers of related 
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responses (Qiow 8 and Qhigh 5 ), wherein each response being raised to a power being a 
positive or negative real number, or zero. 

In preferred embodiments of third aspect of the present invention, the calibrated 
5 functional relationship comprises terms being of the form: Q Iow n1 * Qhigh m1 wherein n1 and 
m2 are real numbers and/or integers, and wherein n1 is positive. In order to for instance 
increase the accuracy of the calibrated functional relationship this relationship may 
preferably comprises terms being of the form: Q !0 w , Qhigh, Qiow 2 , Qwgh 2 and Q !ow / Qhigh, or 
preferably comprises terms of the form: Q| 0W * Qhigh; Qiow 2 *Qhigh; Qiow * Qwgh 2 ; Qiow 2 /Qhigh; 

10 Qiow / Qhigh 2 ; Qiow 2 /Qhigh 2 ; Qiow 2 / Qhigh 2 - 

In particular preferred embodiments according to the third aspect of the present invention 
the calibrated functional relationship is of the form: 

F(fi) = B 1 Q? w + B 2 Q* igh +B 3 QL 2 +B 4 Q $ Mgh 2 +B 5 Qf 0W Q^ igh + B 6 Q* ow 2 Q s high + B 7 Qf ow Q* igh 



15 



, Qiow f p Qiow , t> Vlow t p 

+ B »7^ +B 97^ + B H> c 2+** n 

Vhigh Q 



Clow 



high 



2 

+ B fl 



Vbigh vwgh Vhigh 

It may furthermore be preferred that the calibration model is obtained by exposing the 
master instrument, such as an instrument having carried out the method according to the 
first aspect of the present invention, to a plurality of well-defined objects. 

20 Typically and preferably, well-defined objects are defined in the sense that the physical 
properties of the objects have been established by a chemical process, such as an 
officially recognized reference method for the determination of the requested physical 
properties. 

25 In preferred embodiments and as disclosed in relation to the other aspect of the present 
invention, each of the responses (Q) is preferably either. 

an intensity (I), if necessary corrected with respect to dark current of the 

detectors; 

a transmittance (T) being derived from intensity as a ratio between intensity 
30 resulting from measuring on an object and a reference intensity; 

an absorbance being defined as the negative logarithm to a transmittiance (A=- 
log(T)) such as logarithm base 10, the natural logarithm, or any other logarithmic 
function; 

or 
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a reflectance (R) expressing the reflectance from the surface of the object, the 
reflectance (R) is preferably linearized using the Kubelka-Munk transform 
(K/S=(1-R)/2R). 

In particular preferred embodiments according to the third aspect of the present invention 
where the responses are absorbance, such absorbances are preferably being 
determined by calculating 



A iow = -logio 



psample flOW) -1^ (low) " 

L i air flow)-i djnk aow) _ 



and 



A U8h = -log 



10 



l^Chi^-l^Chigh)] 



1 UChi^-l^Chigh) J 

10 wherein the intensities (l)*are obtained in a measuring region of the slave instrument by: 
exposing an object in the measuring region to low and high X-ray energies and 
detecting with detectors the intensities Umpieflow) and Isampie(high) respectively 
detecting with the detectors the intensities Idark(low) and Idark(high) from said 
detectors when no radiation reaches them; 

15 and 

exposing said detectors to the low and high X-ray energies signals when no 
object is present in the measuring region and detecting l a i r (low) and I a ir(high). 

In a fourth aspect the present invention relates to a method of using a slave instrument 
20 for determining physical quantities. Such physical quantities are for instance: preferably 
the fat content of an object, which object is for instance food or feed, and the quantities 
are preferably determined by use of dual X-ray radiation. In accordance herewith, the 
method comprises preferably: 

scanning substantially all or all of the object by X-ray beams having at least two 
25 energy levels, including a low level and a level being higher relatively thereto, 

detecting the X-ray beams having passed through the object for a plurality of 

areas of the object; 

for each area of the object 

- determining the object's response (Qiow) at the low energy level and the 
30 object's response (Qwgh) at the high energy level, 

- correcting the responses so determined preferably by utilizing the 
correcting method according to the second aspect of the present 
invention, 

and 
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- determining the physical property preferably by utilizing the method according 
the third aspect of the present invention. 

In a sixth aspect of the present invention, the invention relates to a data processing 
5 system for providing a correction for a slave instrument. Such a system utilizes 
preferably sets of responses being based on detecting attenuation and/or reflection 
and/or scatter of electromagnetic radiation, in particular X-ray, in/from a object exposed 
to said electromagnetic radiation in at least two spectral ranges. The set of responses 
comprises preferably one or more pair of related responses (Qiow s and Qhi g h s ) 
10 representing measurements performed with a slave instrument and a set of responses 
comprising one or more pair of related responses (Qiow m and Q h igh m ) representing 
measurements performed with a master instrument. These responses being preferably 
obtained for a plurality of stable objects and 

to each pair of related responses of the slave instrument corresponds a 
15 pair of related responses of the master instrument, 

and to each element in each pair of responses of the slave instrument 
corresponds an element in the corresponding pair of responses of the - 
master instrument. 

20 In accordance with the sixth aspect the data processing system comprises preferably 

means for accessing memory means wherein the responses (Qi ow m and Q h i g h m ) of 
the master instrument and/or the responses (Q Iow s and Q W gh s ) of the slave 
instrument are stored, 

means, such as processor means, for determining based on the sets of 
25 responses a correcting function being a functional relationship between a ratio of 

related responses of the master instrument and a sum of a plurality of terms, 
each term being a product of a correcting coefficient (B|) and powers of related 
responses (CW 5 and Qhigh*) of the slave instrument wherein each response 
being raised to a power being a positive or negative real number, or zero, thereby 

30 determining a first set of correcting coefficients (B 0 ; B 2 ) being multiplied on 

each of the terms, 

means for storing the first set of correction coefficients (B 0 ; B^ B2 ). 



35 



Data processing systems of the sixth aspect of the invention comprise preferably means 
for determining a further correcting function being a functional relationship between a 
high energy response of the slave instrument (Qmgh 5 ) and related high energy response 
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(Qhigh 1 *) of the master instrument, thereby enabling the system to determining a second 
set of correcting coefficients (a; P). 

In accordance with preferred embodiments of the data processing system according to 
the present invention, and as disclosed in connection with the other embodiments of the 
invention, it is preferred that each of the responses (Q) is either 



a transmittance (T) being derived from intensity as a ratio between intensity 

resulting from measuring on an object and reference intensity; 

an absorbance being defined as the negative logarithm to a transmittance (A— 

logCO) such as logarithm base 10, the natural logarithm, or any other logarithmic 

function; 

or 

a reflectance (R) expressing the reflectance from the surface of the object, the 
reflectance (R) is preferably linearized using the Kubelka-Munk transform 
(K/S=(1-R)/2R). 

In particular preferred embodiments of the data processing system according to the 
present invention, in case the responses are absorbances, the absorbances being 
preferably determined by calculating 



wherein the intensities (I) are obtained in a measuring region of the slave instrument by: 
exposing an object in the measuring region to low and high X-ray energies and 
detecting with detectors the intensities Isampie(low) and Isampie(high) respectively 
detecting with the detectors the intensities Idark(low) and ^(high) from said 



an intensity (I), if necessary corrected with respect to dark current of the 



detectors; 





detectors when no radiation reaches them; 



and 



exposing said detectors to the low and high X-ray energies signals when no 
object is present in the measuring region and detecting lairflow) and lair(high). 



In a seventh aspect, the invention relates to a correcting system comprising a slave 
instrument for obtaining responses and a data processing system for correcting the 
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responses, the responses representing measurement performed with the slave 
instruments and the responses being based on detecting by the slave instrument 
attenuation and/or reflection and/or scatter of electromagnetic radiation, in particular X- 
ray, in/from a object exposed to said electromagnetic radiation in at least two spectral 
ranges, the set of responses comprises one or more pair of related responses (Q low s and 
Qhigh s )- In accordance herewith, the correcting system comprises preferably 

processor means for determining the one or more pair of related responses (Q| 0W S 
and Q h igh s ) based on measurement on an object with the slave instrument, 
means comprising processor means adapted to perform a correction of 
responses by utilizing a correcting according to the second aspect of the present 
invention, said processor means comprises 

means for accessing memory means storing a first set of correction 
coefficients (B 0 ; B 2 ) 

processor means for determining the ratio [CW O^] 00 " by the correcting 
function; 

processor means for determining the corrected high energy response 
Qhigh 00 * by the further correcting function; 



and 



processor means for determining the corrected low energy response 
Qiow 00 " by multiplying [Q l0W / by Q high ° 



corr 



Also in this aspect of the invention, and as disclosed in connection with the other aspects 
of the invention, each of the responses (Q) is preferably either: 

an intensity (I), if necessary corrected with respect to dark current of the 

detectors; 

a transmittance (T) being derived from intensity as a ratio between intensity 

resulting from measuring on an object and reference intensity; 

an absorbance being defined as the negative logarithm to a transmittance (A=- 

'og(T)) such as logarithm base 10, the natural logarithm, or any other logarithmic 

function; 

or 

a reflectance (R) expressing the reflectance from the surface of the object, the 
reflectance (R) is preferably linearized using the Kubelka-Munk transform 
(K/S=(1-R)/2R). 
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Preferably, the system according to the seventh aspect comprises storage means 
wherein a set of responses (CW" and Q^O for a set of stable objects measured on a 
master instrument are stored and /or storage means wherein the first set of correction 
coefficients (B 0 ; Bi; B 2 ) and/or the further correcting function is/are stored. 

In an eighth aspect, the present invention relates to a dual X-ray instrument comprising a 
system according to the seventh aspect being adapted to carry out a method according 
to the first aspect of the invention. 



In a ninth aspect, the present invention relates to a set of objects comprising one or more 
stable objects for, or used during, carrying out one or more of the methods according to 
the present invention. Preferably, each of such objects being characterized by being 
composed by at least two different chemical compositions which are substantially stable 
and each stable object is having response, such as absorbance, properties which are 
similar to the response, such as absorbance, properties of an object subjected to the 
method according to the second aspect of the present invention. 

In accordance with the ninth aspect, it is preferred that for each of the stable objects a 
first member of the at least two different chemical compositions is one having X-ray 
response properties, such as absorbance properties, similar to adipose tissue, and a 
second member of the at least two different chemical compositions is one having X-ray 
response, such as absorbance, properties similar to muscle tissue. 

In preferred embodiments of present invention, the set of stable objects comprises 
preferably a plurality of stable objects having varying thickness and/or areal density. The 
plurality of stable objects may advantageously be integrated into a single stepped item. 

Preferably, each object comprised in the set of objects is stable in the sense that the X- 
ray response, such as absorption, properties of the object does not change more than 
0. 1 %, such as no more than 0.01 %, such as no more than.0,001 % within at least 10 
days, such as at least 1 month, preferably at least 1 year. 

Preferably, the number of stable objects comprises in the set of stable objects are at 
least 8, such as at least 12, preferably at least 15, or even at least 20, such as at least 
26. 
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In the following a particular preferred embodiment of the present invention will be 
presented as non-limiting example with reference to the accompanying figures, in which: 

Figure 1 shows a schematic diagram of an instrument incorporating means for 
5 carrying out a method according to the invention. 

Figure 2 shows a perspective view of an embodiment of the instrument in Figure 1 . 
Figure 3 shows an alternative spectrum of X-ray sources simulating two sources by 

use of one source and a combination of two filters. 
Figure 4 shows a typical meat sample in a plastic container. 
10 Figure 5 shows a typical low energy X-ray transmission image of a meat sample as 
shown in figure 4. 

Figure 6 shows a typical high energy X-ray transmission image of the same meat " 
sample. 

Figure 7 is an image illustrating a calculated areal density for each individual pixel. 
1 5 Figure 8 is an image illustrating a calculated fat content for each individual pixel. 

Figure 9 is an image illustrating a calculated "fat map n for a meat sample of 36 % fat. 
Figure 10 shows an embodiment of an integrated standardization object. 
Figure 1 1 shows an embodiment of individual steps of an integrated standardization 
object. 

20 Figure 12 shows a second embodiment of an integrated standardization object. 
Figure 1 3 shows a further embodiment of a number of steps and sections of an 

integrated standardization object. 
Figure 14 shows a flow diagram illustrating the measuring process with a master 

instrument 

25 Figure 1 5 shows a flow diagram illustrating the standardisation process. 

Figure 16 shows a flow diagram illustrating the new measurement process using a 
standardized slave. 

Figure 1 7 shows a plot of ratios of non-standardised absorbances obtained with four 
non-standardised instruments against ratios of absorbances obtained with 
30 a master instrument, 

Figure 1 8 shows a plot of ratios of standardised absorbances obtained with four 
standardised instruments against ratios of absorbances obtained with a 
master instrument, 
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DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT OF THE PRESENT 
INVENTION 

The present invention may be applied for instruments of the kind measuring 
absorbances, in a medium, of an electromagnetic radiation in at least two spectral 
ranges, which instrument is calibrated by use of a multivariate regression method. 
Such instruments may be spectral IR-instruments. However, the present invention has 
proved to be specifically useful for Dual X-ray analysis. 

Description of the equipment used 

The following description discloses as an example a preferred embodiment of an 
instrument for which the present invention is specifically intended. The instrument uses 
two X-ray sources. The instrument is designed for being installed in relation to a 
production iine in a slaughterhouse. Figure 1 shows a schematic diagram of an 
embodiment of a measurement system for a determination of the fat content in meat. 
Figure 2 is a perspective illustration of the presently preferred X-ray instrument. Figure 2 
shows only the active operating portions of the X-ray equipment. For purpose of clarity, 
all protective shielding or screening and all casings are deleted from the drawing. The - 
equipment comprises or is located in close relation to a conveyor 10. Two X-ray sources 
12, 14 are arranged above the conveyor 10. From the two sources 12, 14 X-ray beams 
1 6, 1 8 are directed towards detectors 22, 24 arranged below the conveyor. The 
conveyor may be split into two separate conveyors spaced to allow free pass of the X- 
rays and to leave an open space for location of detectors 22, 24. Alternatively the 
conveyor belt should be made from a material showing a low absorbance of X-rays, e.g. 
polyurethane or polypropylene. The food or feed to be measured is arranged in an open 
container or box 20, preferably also composed by a material showing low absorbance of 
X-rays. Alternatively the medium such as a food product might be arranged directly on a 
conveyor belt. In a further alternative arrangement the X-ray sources could be located 
below the conveyor and the detectors above the conveyor. 

The presently preferred equipment used in the present example consists of two constant 
potential X-ray sources 12, 14, one at low energy (e.g. 62 kV/5.5 mA) and another at 
high energy (e.g. 120 kV/3.0 mA), both with an appropriate filtration (e.g. using 0.25 and 
1.75 mm of copper, respectively) narrowing the spectral range of the radiation emitted 
from the polychromatic sources. The two sources are spatially separated to avoid 
interference between them, i.e. to avoid that radiation from one source is detected as if it 
originated from the other. The radiation from either source is collimated by a lead 
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collimator. In this way two fan-shaped beams of X-rays 16, 18 are directed through 
container 20 comprising a sample or batch of the food or feed product towards detectors 
22, 24, e.g. Hamamatsu C 7390. Alternatively the meat lumps may be arranged loosely 
on a conveyor band. 

Further, the two separate sources may be replaced by a combination of one source and 
two filters emitting a low energy and a high energy beam. The resulting source spectra 
are shown in Figure 3. However the preferred embodiment applies two separate sources 
12, 14 driven by separate power supplies 13, 37. 

Each of the two X-ray sources 12, 14 is associated with an array of detectors 22, 24 
covered with a scintillating layer converting the transmitted radiation into visible light that 
can be measured by the detectors 22, 24. The scintillating layer may consist of e.g. 
cadmium telluride, mercury iodide, cesium iodide (Csl), gadolinium oxysulphide 
(Gd 2 0 2 S), or yttrium oxysulphide (Y 2 0 2 S), and/or CdW0 4 , preferably doped in order to 
reduce the after-glow effect. The pixels used in the presently preferred embodiment have 
the dimensions 1.6x1.3 mm 2 and are arranged as an array of 384 pixels with a pitch of - 
1 .6 mm. These dimensions are only stated as an example. Other dimensions may be 
applied. The pixels convert the amount of transmitted light into analogue signals that are 
passed through cables 27, 28 to an analogue-to-digital converter 34, which is connected 
through cable 35 to a computing means 38 capable of performing the successive 
calculations. 

A monitor 42 may be connected through cable 40 to the computing means to show ' 
results or details of the operation. The computing means 38 may include means for 
controlling the supply of power through means 36, 37, 26 and 25, 13, 15 to the X ray 
sources 12, 14. The monitor 42 and the computing means 38 may comprise a Personal 
Computer, preferably including at least one Pentium processor and/or a number of digital 
signal processors. 

The operational speed of the conveyor is preferably substantially constant. The items, 
motor 30, control box 33, and cables 32, 39, shown by phantom lines in Figure 1, 
indicate that the operation of the conveyor optionally may be controlled by the computing 
means 38. The conveyor may include position measuring means, e.g. an encoder 
installed on a conveyor driving shaft. Alternative means may be a laser or radar detection 
or marks on the conveyor belt. It is essential to the present method that the data 
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representing the two X-ray images can be synchronised. Such synchronisation may 
however be obtained in many ways, including mathematical post-processing of the 
images. 

Operation of the instrument during a normal measurement 

A container 20, e.g. as shown in Figure 4, comprising e.g. meat trimmings from a cutting 
section of the slaughterhouse, is received on the conveyor 10. The container is moved 
with a fairly constant speed of e.g. about 5 - 100 cm per second, such as 10-50 cm, e.g. 
30 cm per second past the fan shaped beams 16, 1 8 and the arrays of detectors 22, 24 
in a controlled manner in order to generate two "images" of the absorbances in the 
sample or batch, one at a low X-ray energy (shown in Figure 5) and another at a high 
energy (shown in Figure 6). All data representing the two images are stored in the 
computer 38. 

Treatment of the collected data from a master instrument: 

Fig. 14 represents a flow chart illustrating the measurement and data treatment in a 
master instrument. As stated above, data representing two X-ray images of (Figures 5, 6) 
of each container (Figure 4), comprising a batch of food or feed e.g. meat, are obtained. 
The signals at the pixels are l low and at low and high X-ray energies, respectively, 
(1 10, 1 12 in Figure 14). Furthermore, the so-called "dark signals" (i.e. the signal from the 
detectors when no radiation reaches them), l da *(low) and l dartc (high), and the "air signals" 
(i.e. the signal from the detectors when no sample is present in the measuring region), 
lair(low) and l air (high), are collected for each pixel at both X-ray energies (102 in Figure 
14). Preferably these data are collected repetitively in the intervals between the 
passage/passing of meat containers, i.e. the dark signals and air signals are measured 
repetitively, e.g. at regular intervals during a day to adjust for any drift of instrument 
performance. 

Now referring to 1 14 in Figure 14, these signals are transformed into absorbance units 
by using the following formulas: 

a _ w P—jae flow) -Inflow)] 

A, ow - -log 10 p — — — .(ia) . 

L lair GOW)- 1 ^ (l0W) J 
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A -J _ _ I sample 

m ' ~ l0 n u^-i^chigh) J <ib) 

If the instrument were a calibrated master instrument the arear density (Figure 7) and fat 
content (Figure 8, 9) of the measured medium could now be calculated from these 
absorbances by use of the fat calibration as described in WO 01/29557 and shown 
schematically in box 116 in Figure 14. 

Description of the method and means according to the present invention: 
A.: Description of the standardisation objects. 

Figures 10-13 show various preferred embodiments of integrated standardization 
objects. The present method requires a plurality of stable objects, also called artificial 
samples due to the choice of material, which must be of a kind maintaining substantially 
constant absorption properties for a long period, such as several years, contrary to the 
medium, such as meat, for which this method is specifically - but not exclusively - 
intended.. Preferably, the stable objects include a number of combinations of two X-ray 
absorbers, having absorption characteristics similar to the measured object, which in the 
present example is meat. Accordingly, in the present example, the two absorbers must 
behave like adipose and muscle tissue. In the presently preferred embodiment blocks of 
polymethyl methacrylate 81 (similar to adipose tissue) and "plastic water" 82 (a polymer, 
with absorption characteristics similar to water, available from CIRS, Inc.) were used. 
Other alternatives, such as other purpose-made polymers made to resemble adipose 
and muscle tissue, could also be used. Liquid water (or ice) is another substance highly 
resembling muscle tissue. 

In the present example twenty-six samples consisting of various combinations of these 
materials 81 82 were generated (areal densities ranging from approx. 1 to 22 g/cm 2 ) and 
used in the standardisation procedure. 

In order to make the standardization process easy to perform it is preferred that the 
whole number of different stable objects are integrated into a single item, e.g. a stepped 
item, e.g. like a staircase 50 as shown in Figure 10 and 12. As shown in Figures 11 and 
1 3 each step 51 , 52, 70, 71 ... .79 may comprise one or two or an other plurality of 
different layers 81, 82 of e.g. the polymethyl methacrylate and "plastic water". Further 
each step 51, 52, 70 may comprise several sections 53, 54, 55, ...59, each having a 
specific combination of layers 81, 82 of various thickness. Preferably, the different 
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combinations are chosen to provide a range of absorptions and properties, e.g. fat % 
similar to the absorptions and properties, e.g. fat%, in the media for which the instrument 
is intended. 



5 Preferably, each section has the same size in the travel direction of the conveyor. 
Preferably, they also have the same size in a direction perpendicular to the travel 
direction and parallel to the conveyor belt. In other word the projection area on the 
generally horizontal plane of the conveyor belt of each section is the same in a preferred 
embodiment. 

10 

Each standardisation object is measured on a master instrument or standardized slave 
instrument. Preferably, each standardization object is marked by an identification code, 
such as a number. The measured absorbances and preferably the identification code are 
stored in a memory, which is accessible from a slave instrument, using the 
15 standardization object. 

B standardization of a slave instrument 

The standardisation procedure requires a set of absorbances, A\w m and Awgh m , from a 
master instrument, and a set of corresponding absorbances, A| OW s and Ahi 9 h s , from a slave 
20 instrument. These absorbances must originate from the same samples, e.g. a number of 
artificial standardisation samples, preferably being a standardization object as described 
above. 

A presently preferred method according to the invention is shown schematically in the 
25 diagram Figure 15. In step 202 l dar k(low) and l d ark(high),and Uflow) and l a ir(high), are 
(collected for each pixel at both X-ray energies. Preferably these data are collected 
repetitively in the intervals between the passage/passing of meat containers, i.e. the dark 
signals and air signals are measured repetitively, e.g. at regular intervals during a day to 
adjust for any drift of instrument performance. 

30 

In step 204 a standardization object 50 is arranged on the conveyor 10 (shown in Figure 
1 and 2). The object is passed through the instrument in the same manner as an ordinary 
object to be measured. The steps 200, 202, 204, 206, 208, 210, 212, are the same as 
mentioned before regarding Figure 14. 
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In the following step 214 average values /W, Aw g h s are calculated for each section 53, 
54, 55, ...59, of the standardization object 50 measured on the slave instrument. The 
values are calculated using formula (1a) and (1b). The high and low energy values are 
coordinated, i.e. the values belonging to the same section (53, 54, 55, ...59, ) are 
5 matched. The corresponding values A low m , A h i gh m , measured on a master instrument are 
read from a memory. 

The following final steps listed in box 216 are then required to obtain an acceptable 
standardisation of the instruments: 

10 

1 . The coefficients (a and p) in the equation describing the relationship between A Wgh m 
and Ahigh s are determined by univariate linear regression: 

This step may in some cases be superfluous as a may be almost equal to 1 and p close 
15 to zero. 

2. A functional relationship, /(A| 0 w s ,Ahigh s ), between the A, ow m /Ah Igh m ratio and various 
combinations of A, ow s and Aw gh s is established. The following linear combination is 
preferred: 

^ = B,Af 0W + B 2 + B 3 A* J + B 4 A s high 2 + B 5 A^A^ + B 6 A^A^ + B 7 A^ A^ 2 

^high 

20 s 2 s r s I 2 

+ B 8 -^ + B 9 4|^ + B 10 ^V + B I1 i±js-\ +B 0 =/(A? OT ,A^) 
Ahjgh L A ^J 

The coefficients, B 0 to Bn , are be determined e.g. by using Partial Least Squares 
(PLS) regression. Alternatively other multivarialte regression methods may be used. 

3. The coefficients, a, p, and B 0 to Bn are saved in the instrument. They are the 
25 instrument specific standardisation (correction) coefficients. 

C. Measuring an unknown meat sample on the standardized slave instrument: 

A presently preferred procedure for measuring a new meat sample on the slave 
30 instrument is illustrated in Figure 1 6. A container (Figure 4) is arranged on the conveyor 
1 0. During the passage through the instrument the two X-ray beams pass through the 
meat, container bottom, conveyor belt arid reach the detectors 22, 24, generating signals 
representing the two images shown in Figures 5, 6. Data representing the intensities li™, 
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l hlgh are stored temporarily. I da rk and I* are measured regularly in step 312. Absorbances 
Atow S , A hlgh s are calculated by use of formula (1a) and (1b). For each pixel -or 
alternatively a small group of neighboring pixels - two matching values Ajow 8 , Ahi 9 h s are 
coordinated. So far the data treatment is the same as shown in Figure 14. 
5 According to the present invention the following treatment is applied to the raw 
measurements, A,™ 8 and A W3 h s , from the slave instrument: 

1 . The corrected high energy absorption, Awgh 00 "; is calculated: 

A h!gh corf =a.A h igh S +P 
10 by using a and p determined as disclosed above in section C. 



2. The corrected absorption ratio, [A to JA h jgh] corr , is calculated: 



[ 



7 ~ /(Alawy ^high) 



k high 



15 where/ is the function defined above in section C. 



3. The corrected low energy absorption, Aiow 00 ", is calculated: 



A corT — 



Mow 



A coir 
^high 



20 These corrected absorbances, A,™ 00 "" and A^ 00 ", are hereafter used for predictive 
purposes, e.g. prediction of the fat content of a meat sample using a calibration model 
generated on the master instrument (or any other instrument in the population of 
standardised instruments). This is indicated in the box 316, which is identical to box 116 
in Figure 14. 

25 

The example presented below shows how the method may work in practice. 

Meat samples used for exemplifying the method according to the invention: 
86 samples consisting of minced pork, turkey, and beef meat were prepared. These 
30 samples were frozen in blocks of varying heights (from 1 0 to 200 mm, corresponding to 
areal densities from 1 to 21 g/cm 2 ). The fat content of these samples was determined 
using the fat reference method (SBR, Schmid-Bondzynski-Ratzlaff). It ranged from 2 to 
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73 %. 44 of these samples were used for building calibration models while the remaining 
42 samples were used for independent testing of the standardisation method. 



Description of the measurements 

5 The 26 artificial samples and the 86 meat samples were measured on equipment as 
described above. According to previous experience it is known that particularly voltage 
changes in the low X-ray sources cause serious problems with the transferability of the 
calibrations. Therefore, the experiment was carried out five times with different voltage 
settings of the two energy sources, on different instruments. The following settings were 

io used: 



Instrument 1 (master): 


E|0W 




62 


kV, 


Ehlgh = 


120 


kV 


Instrument 2 (slave): 






58 


kV, 


Ehfgh = 


100 


kV 


Instrument 3 (slave): 


E|ow 




58 


kV, 


Ehlgh = 


110 


kV 


Instrument 4 (slave): 


Elow 




66 


kV, 


Ehlgh = 


110 


kV 


Instrument 5 (slave): 


Elow 




66 


kV, 


Ehlgh = 


100 


kV 



Use of the method 

PLS calibration models for fat and areal density based oh 44 meat samples werie 
20 calculated on data from Instrument 1 . This instrument is therefore regarded as the 

master instrument against which the so-called slave instruments (Instruments 2 to 5) will 
be evaluated. The calibration models consist of linear combinations of various ratios and 
products of A !ow and A h i g h, such that the fat % in a specific specific area or point of an 
object, above a specific single detection element, or alternatively a neighborhood of 
25 adjacent detection elements, such as four or nine, may be calculated as: 

Fat%= b 0 + bf *A l0W + b 2 *A hj g h + b 3 *A, 0W 2 + b 4 * A hlgh 2 f b 6 * A, ow /A hlgh /+ 

. b 6 * A,™* A h!gh + br * A low 2 * A high + b 8 * A, ow * A hIgh 2 + b 9 * A l0 w* A high ^ + 
bio* Ajow^Ahigh 4 +b 11 *A,ow 2 /A Wgh +b 12 *A low /A h!gh 2 + b 13 * Aow 2 /^ 2 
30 + b 14 * A IOW 3 / A hlgh 2 + b 15 * A l0W 4 / A hjgh 2 + b 16 * 1 /Ah igh 4 + b 17 * A, 0W 4 / A high 3 + 
bis * Atow 3 / Ahigh 4 + b 19 * Aiow 4 / Ahigh 4 + + bp * Aiow n1 / Ahigh" 11 ., 

wherein b 0 , b 1t ... b p (some of them may be zero) are the calibration coefficients 
determined through multivariate calibration. Typically, the series is truncated so as to 
35 contain only term up to power two, such as power three, preferably power four of A: 
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It is contemplated that the expected accuracy of a calibrated instrument against the 
reference method (expressed as the Root Mean Square Error of Prediction, RMSEP) is 
better than 1 % for the fat determination and approx 0.1 gfcm 2 for the areal density 
determination. 

5 

The A !cw /Ahigh ratio is the parameter that is most sensitive to whether the instrument is 
standardised or not. This is of major importance, since this ratio is one of the major 
contributors to the calibration model for fat. In Figure 17 this ratio for the 86 meat 
samples measured on the master instrument is plotted against the non-standardised 

10 ratios for the slave instruments. If no standardisation problems existed, the points should 
be close to the straight line indicated in the figure. This is, however, far from true: 
especially Instrument 5 shows large errors. A simple slope and intercept correction of 
this error will not help the problem completely, as the points, apart from showing a large 
systematic error, are also scattered along an imaginary line passing through the cloud of 

15 points. This is the reason why eleven or more terms are required in f(A| 0 w s ,Ah]gh s )- 

After calculation of the instrument dependent standardisation (correction) coefficients; ar 
p, and B 0 to Bn, for Instrument 2 to 5 from the data obtained on the 26 artificial samples, 
the absorbances of the 86 meat samples were corrected using these coefficients. This 
20 resulted in a set of corrected ratios, [Aiow/Awgh] 00 ", that are plotted in Figure 18. The ratios 
are now very close to the line, indicating that the ratio is independent of the instrument 
from which it originated, which means that calibration models can be transferred between 
instruments. 

25 The result of various attempts to use the calibration models based on data from the 

master (Instrument 1) are presented in Table A. Firstly, the calibration models for fat and 
areal density were applied to the raw, non-standardised data from the 42 independent 
test meat samples (corresponding to the data presented in Figure 17). As is evident from 
column two and four of Table A, this results in very large and unacceptable prediction 

3d error. When standardisation is applied (corresponding to the data presented in Figure 
1 8), on the other hand, the prediction errors for the slaves (Instruments 2 to 5) cannot be 
distinguished from the prediction error for the master (Instrument 1). This can be seen in 
columns three and five of Table A. 



35 



# + 
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Non- Norv 

standardised Standardised standardised Standardised 

fet fet area) density area) density 

prediction error prediction error predicfion error predicfon error 

Instrument (%) (%) (gfcrrf) (gfcnf) 



1 master 


0.66 


0.66 


0.10 


0.10 


2 slave 


4.56 


1.06 


0.43 


0.12 


3 stave 


20.69 


0.90 


0.32 


012 


4 slave 


38.08 


0.72 


019 


0.11 


5 slave 


54.07 


1.15 


O09 


0.11 



(RMSEP) 
Comments 

5 The example given above thus shows the advantages of using the described method for 
standardising an X-ray instrument The major advantage lies typically in the fact that it is 
possible to use stable artificial samples with well-defined absorption characteristics for 
obtaining standardised fat and areal density predictions instead of having to perform a 
labour intensive and expensive calibration procedure for every new instrument or every 

10 time an instrument is modified, e.g. by replacing a defect X-ray source or detector. 

While a single particular embodiment of the invention has been mentioned, it will be 
understood, of course, that the invention is not limited thereto since many modifications 
may be made. It is contemplated that it will be useful for other kinds of electromagnetic 
15 radiation as well , e.g. Infrared light. It is, therefore, contemplated to cover by the 

appended claims any such modifications as fall within the true spirit and scope of the 
invention. 



m * 
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CLAIMS 

1. A method of providing a correcting for a slave instrument, of the kind measuring 
properties of an object by exposing the object to electromagnetic radiation, in particular 
5 X-rays, in at least two spectral ranges and obtaining one or more object responses 
thereto, the responses being based on detecting attenuation and/or reflection and/or 
scatter of the electromagnetic radiation in/from the object by use of one or more 
detectors and are obtained in a form where they express properties either directly or via 
a transformation, said method of correcting comprising 
10 - obtaining, for a plurality of stable objects, a set of responses comprising one or 
more pair of related responses (Q| OW s and Q h i 9 h s ) representing measurements in 
the at least two spectral ranges performed with the slave instrument and a set of * 
responses, comprising one or more pair of related responses (Qkw m and QhigO 
representing measurements in the at least two spectral ranges performed with a 
15 master instrument, 

to each pair of related responses (Q loW s and Q Wgh s ) of the slave instrument 
corresponds a pair of related responses (Q low m and Qhi 9 h m ) of the master - 
instrument, 

and to each element in each pair of responses (Qjow s and Q W gh s ) of the 
20 slave instrument corresponds an element in the corresponding pair of 
responses (Qiow m and Q W gh m ) of the master instrument; 
- determining based on the sets of responses a correcting function being a 
functional relationship between a ratio of related responses of the master 
instrument and a sum of a plurality of terms, each term being a product of a 
25 correcting coefficient (B f ) and powers of related responses (Q low s and Qwgh^ of 
the slave instrument, wherein each response being raised to a power being a 
positive or negative real number, or zero, thereby determining a first set of 
correcting coefficients (B 0 ; Bi; B 2 ) being multiplied on each of the terms; 

and 

30 - storing the first set of correcting coefficients (B 0 ; B^ B 2 ) in memory means 

included in or adapted for communication with data processing means included in 
or adapted for communication with the slave instrument 



35 
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2. A method according to claim 1, wherein, initially, at a manufactures site 

measuring the plurality of stable objects on a master instrument, thereby 
obtaining the set of responses representing measurements performed with the 
master instrument (Qiow m and QhigO, 
5 - storing the set of responses (Qiow m and QhigO as a set of constant values in 

memory means, which is accessible from a slave instrument, when measuring the 
corresponding stable objects on a slave instrument in order to carry out a method 
of providing a correcting according to claim 1. 

10 3. A method according to claim 2, wherein the set of responses measured by the master 
instrument is stored in memory means included in or adapted for communication with 
data processing means included in or adapted for communication with the slave 
instrument 

15 4: A miethod according to claim 1 , 2 or 3, wherein the determination of the correcting 
function being based on a regression method. 

5. A method according to claim 4, wherein the regression method is selected from the 
group consisting of principal component regression, multiple linear regression, partial 

20 least squares regression, and artificial neural networks. 

6. A method according to any of the preceding claims, wherein the correcting function 
comprising a plurality of terms of the following form Qiow n1 * Qh>gh m1 wherein n1 and rn2 
are real numbers and/or integers, and n1 is positive. 

25 

7. A method according to claim 6, wherein the correcting function comprising at least 
three of the following terms: Q ]aw , Q h i gh , Q IOW 2 , Q h!gh 2 and / Q hIgh . 

8. A method according to claim 6, wherein the correcting function comprising at least 
30 three of the following terms: Q,^ * Q^; Q low 2 *Q high ; Q IOW * Q hlgh 2 ; Q !ow 2 /Q h ign; Qiow / Qwgh 2 ; 

Qtow 2 /Qhigh 2 ; Q|ow 2 /Qhlgh 2 - 
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9. A method according to any of the preceding claims, wherein the correcting function is 
of the form: 

^ = B,QL +B 3 QL 2 +B 4 Q^ 2 +B 5 Qr 0W Q 5 hieh +B 6 Q? 0W 2 Q 8 W8h +B 7 Q? 0W QL 8h 5 

r -12 



Vhigh 



10. A method according to any of the preceding claims, further comprising determining 
based on the sets of responses a further correcting function being a functional 
relationship between responses of the slave instrument (Qiow 3 or Qhigh 3 ) and related 
responses (Q !o w m or Q W9 h m ) of the master instrument, thereby determining a second set 

10 of correcting coefficients (a; p). 

1 1 . A method according to claim 1 0, wherein the further correcting function being a 
functional relationship between a high energy response of the slave instrument (Qwgh s ) 
and the related high energy response (Qhi gh m ) of the master instrument. 

15 

12. A method according to claim 1 1 , wherein the further correcting function is determined 
by use of univariate linear regression. 

13. A method according to claim 12, wherein the further correcting function being of the 
20 form Q W gh m = a-Qhigh 5 + p. 

14. A method according to any of the preceding claims, wherein the set of responses for 
the master instrument and the set of responses for the slave instrument each comprise 
one pair of related responses for each stable object comprised in the plurality of stable 

25 objects. 

15. A method according to any of the preceding claims, wherein the related responses 
are obtained based on measuring on objects being conveyed. 



30 



16. A method according to any of the preceding claims, wherein each of the responses 
(Q) is an intensity (I), if necessary corrected with respect to dark current of the detectors. 
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17. A method according to any of the claims 1-15, wherein each of the responses is a 
transmittance (T) being derived from intensity as a ratio between intensity resulting from 
measuring on an object and reference intensity. 

5 1 8. A method according to any of the claims 1-1 5, wherein each of responses is an 
absorbance being defined as the negative logarithm to a transmittance (A=-log(T)) such 
as logarithm base 10, the natural logarithm, or any other logarithmic function. 

* 

19. A method according to any of the claims 1-15, wherein the responses for both the 
10 master and the slave instruments are absorbances being determined by calculating 



A low =~log 



10 



sample 

I^flow)-!^™) J 



and 



a 11 sample 

(high) -I^ (high)] 
- l0 g"{ ^(hi^-l^Chigh) J 

wherein the intensities (I) are obtained in a measuring region of the master instrument 

respectively the slave instrument by: 
15 - exposing the object in the measuring region to low and high X-ray energies and 
detecting with detectors the intensities Isampie(low) and l 5 am P ie(high) respectively 
detecting the intensities l da rk(low) and Idark(high) from said detectors when no 
radiation reaches them; 

and 

20 - exposing said detectors to the low and high X-ray energies signals when no 
object is present in the measuring region and detecting Inflow) and l a ir(high). 

20. A method according to any of the claims 1-15, wherein each of the responses is a 
reflectance (R) expressing the reflectance from the surface of the object. 

25 

21. A method according to claim 19, wherein the reflectance (R) is linearized, preferably 
using the Kubelka-Munk transform (K/S=(1-R)/2R). 

22. A method of correcting responses representing measurements performed with a 
30 slave instrument, said method comprising for an object 

determining based on measurements with the slave instrument a pair of related 
responses (Q low s and Q h ig h s ), 

determining the ratio [Q lovv / Qwgh] 00 "" by a correcting function being a functional •; 
relationship between a ratio of related responses of the master instrument and a 



10 



• • 1 
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sum of a plurality of terms, each term being a product of a correcting coefficient 
(B,) and powers of related responses (Q low s and Q W9h 8 ) of the slave instrument 
wherein each response being raised to a power being a positive or negative real 
number, or zero, 

providing Q m ™ n either by assuming that Qhigh 00 "' is substantially equal to or 
by use of a further correcting function correlating Qhigh 001 * with Qwgh 8 ; 

and 

calculating CW^ as Qhigh 00 " * [CW QhWiT" I 
thereby providing a set of corrected responses. 

23 A method according to claim 22, wherein the further correcting function being of the 
form Qhigh 00 " = ct-Qhigh 8 + P- 

24. A method according to claim 22 or 23, wherein the correcting function comprises 
15 terms of the following form Q low n1 * Q h igh m \ wherein n1 and m2 are real numbers and/or 

integers, and wherein n1 is positive. 

25. A method according to any of the claims 22-24, wherein the correcting function 
comprises at least three of the following terms: Qlow , Qwgh, Qiow 2 , Qhigh 2 and Q )ow / Q W gh • 

20 

26. A method according to any of the claims 22-25, wherein the correcting function 
comprising at least three of the following terms: Qiow* Qhigh; Qiow 2 *Qhign; Qtow * Qhigh 2 ; Qiow 

/Qhighl Qlow / Qhigh 2 *. Qlow 2 /Qhigh 2 ! Qlow 2 / Qhigh 2 

25 27. A method according to any of the claims 22-26, wherein the correcting function is of 
the form: 



Qhigfr > 



= B,Q' ow +B 2 Qu gh +B 3 Q? ow 2 +B 4 QT m 2 +B S QUQ^ +B 6 Q? W 2 Q^ +B 7 QU<W 



+ B 8 ^ + B 9 ^+B J0 -^- + B n 
^Jhigh Vhigh Mbigh 



Q 



low 



Qhigh J 



2 



+ B„ 
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28. A method according to any of the claims 22-27, wherein each of the responses (Q) i 
an intensity (I), if necessary corrected with respect to dark current of the detectors. 
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29. A method according to any of the claims 22-27, wherein each of the responses is a 
transmittance (T) being derived from intensity as a ratio between intensity resulting from 
measuring on an object and a reference intensity. 

5 30. A method according to any of the claims 22-27, wherein each of responses js an 
absorbance being defined as the negative logarithm to a transmittance (A=-log(T)) such 
as logarithm base 10, the natural logarithm, or any other logarithmic function. 

31. A method according to any of the claims 22-27, wherein the responses are 
10 absorbances being determined by calculating 

"Inflow) -JdgkO 0 ^ 



A low =-log 



10 



and 



(low) -I^ (low) J 

A i sample 

(high) -1^ (high) 
A high - 1Ogl0 L I air (Wgh)-I^(hig^ J 
wherein the intensities (I) are obtained in a measuring region of the slave instrument by: 
exposing an object in the measuring region to low and high X-ray energies and _ 
15 detecting with detectors the intensities Isampie(low) and Umpiefhigh) respectively 

detecting with the detectors the intensities l da rk(low) and l d ark(high) from said 
detectors when no radiation reaches them; 

and 

exposing said detectors to the low and high X-ray energies signals when no 
20 object is present in the measuring region and detecting Wow) and UChigh). 

32. A method according to any of the claims 22-27, wherein each of the responses is a 
reflectance (R) expressing the reflectance from the surface of the object. 

25 33. A method according to claim 32, wherein the reflectance (R) is linearized, preferably 
using the Kubelka-Munk transform (K/S=(1 -R)/2R). 



30 



34. A method according to any of the preceding claims, wherein the correcting function 
arid the further correcting function being determined by the method according to any of 
the claims 1-21. 
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35. A method of determining a physical quantity for an object by a slave instrument, the 
method comprising 

determining for the object corrected high and low energy responses (Q h ig h corr and 
Qiow 0011 ) by utilizing the method according to any of the claims 22-34, 
determining the physical quantity by applying on said corrected responses a 
calibrated functional relationship between Qhigh 00 " and Ok,* 00 " and a physical 
quantity. 

36. A method according to claim 35, wherein the calibrated functional relationship being 
a functional relationship between a physical quantity and a sum of a plurality of terms, 
each term being a product of a calibration coefficient (B,) and powers of related 
responses (Q Iow s and Q Wg h S ), wherein each response being raised to a power being a 
positive or negative real number, or zero. 

37. A method according to claim 36, wherein the calibrated functional relationship 
comprises terms being of the form: Qu^' * Q hIgh m1 wherein n1 and m2 are real numbers 
and/or integers, and wherein n1 is positive, such as comprises terms being of the form: 
Qiow , Qhigh, Qtow 2 , Qh.g h 2 and Q lw / Qwg h , preferably comprises terms of the form: CW 

Qhigh; Qlow 2 *Qh)gh; Qiow * Qhigh 2 ! Qlow 2 /Qhlgh; Qiow / Qhigh 2 ! Qiow 2 /Qhigh 2 ; Qiow 2 / Qhigh 2 - 



38. A method according to claim 37, wherein the calibrated functional relationship is of 
the form: 

F02) = B 1 Q, S OW +B 2 Q S ^ +B 4 QL gh 2 +B 5 QLQ« eh ♦B^V** + B 7 Q tow Q 

+ B 8 -— + B 9 -— + B 10 2 +tJ„ +D o 

Qwgh Qwgh QLgh L^^J 

39. A method according to any of the claims 35-38, wherein the calibration model is 
obtained by exposing the master instrument to a plurality of well-defined objects. 

40. A method according to claim 39, wherein the well defined objects are defined in the 
sense that the physical properties of the object have been established by a chemical 
process, such as an officially recognized reference method for the determination of the 
requested physical property. 
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41. A method according to any of the claims 35-40, wherein each of the responses (Q) is 
either 

an intensity (I), if necessary corrected with respect to dark current of the 

detectors; 

5 a transmittance (T) being derived from intensity as a ratio between intensity 

resulting from measuring on an object and a reference intensity; 
an absorbance being defined as the negative logarithm to a transmittance (A=- 
log(T)) such as logarithm base 10, the natural logarithm, or any other logarithmic 
function; 

10 or 

a reflectance (R) expressing the reflectance from the surface of the object, the 
reflectance (R) is preferably linearized using the Kubelka-Munk transform 
(K/S=(1-R)/2R). 

15 42. A method according to claim 41 , wherein; in case the responses are absorbances, 
the absorbances being determined by calculating 



. r i sa ro p1e0ow)-I dark qOW)] 

AW= -'° B 4 U0ow)- Id-0 o W ) J " ,d 



10 



I sample 

(high) -I*,* (high)] 
L I* (high) (high) J 



wherein the intensities (I) are obtained in a measuring region of the slave instrument by: 
20 - exposing an object in the measuring region to low and high X-ray energies and 
detecting with detectors the intensities l S am P ie(low) and l samp ie(high) respectively 
detecting with the detectors the intensities l<jaiic(low) and ldark(high) from said 
detectors when no radiation reaches them; 

and 

25 - exposing said detectors to the low and high X-ray energies signals when no 
object is present in the measuring region and detecting Utfow) and l a ir(high). 

43; A method of using a slave instrument for determining physical quantities, such as the 
fat content, of an object, such as food or feed, by use of dual X-ray radiation, said 
30 method comprising: 

scanning substantially all or all of the object by X-ray beams having at least two 
energy levels, including a low level and a level being higher relatively thereto, 
detecting the X-ray beams having passed through the object for a plurality of 
areas of the object; 
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for each area of the object 

- determining the object's response (Q low ) at the low energy level and the 
object's response (Qw g h) at the high energy level, 

- correcting the responses so determined by utilising the correcting 
method according to any of the claims 22-34, 

and 

- determining the physical property by utilizing the method according to claims 35- 
42. 

44. A data processing system for providing a correction for a slave instrument, said 
system utilizes sets of responses being based on detecting attenuation and/or reflection 
and/or scatter of electromagnetic radiation, in particular X-ray, in/from a object exposed 
to said electromagnetic radiation in at least two spectral ranges, the set of responses 
comprises one or more pair of related responses (Qiow s and Qnigh®) representing 
measurements performed with a slave instrument and a set of responses comprising one 
or more pair of related responses (Q, ow m and Q h i 9 h m ) representing measurements 
performed with a master instrument, said responses being obtained for a plurality of 
stable objects and 

to each pair of related responses of the slave instrument corresponds a 
pair of related responses of the master instrument, 
and to each element in each pair of responses of the slave instrument 
corresponds an element in the corresponding pair of responses of the 
master instrument; 
said data processing system comprising 

means for accessing memory means wherein the responses (Qiow m and Q h jgh m ) of 
the master instrument and/or the responses (Qiow s and Q h]gh *) of the slave . 
instrument are stored, 

means, such as processor means, for determining based on the sets of 
responses a correcting function being a functional relationship between a ratio of 
related responses of the master instrument and a sum of a plurality of terms, 
each term being a product of a correcting coefficient (Bi) and powers of related 
responses (Q Iow s and Q h i 9 h s ) of the slave instrument wherein each response 
being raised to a power being a positive or negative real number, or zero, thereby 

determining a first set of correcting coefficients (B 0 ; 6^82 ) being multiplied on 

each of the terms, 

means for storing the first set of correction coefficients (B 0 ; B t ; B 2 .„...). 
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45. A data processing system according to claim 44, further comprising means for 
determining a further correcting function being a functional relationship between a high 
energy response of the slave instrument (Qwgh 5 ) and related high energy response 
(Qwgh") of the master instrument, thereby determining a second set of correcting 
coefficients (a; p). 

46. A data processing system according to claim 44 or 45, wherein each of the 
responses (Q) is either 



detectors; 

a transmittance (T) being derived from intensity as a ratio between intensity 

resulting from measuring on an object and reference intensity; 

an absorbance being defined as the negative logarithm to a transmittance (A— 

Iog(T)) such as logarithm base 10, the natural logarithm, or any other logarithmic 

function; 

or 

a reflectance (R) expressing the reflectance from the surface of the object, the 
reflectance (R) is preferably linearized using the Kubelka-Munk transform 
(K/S=(1-R)/2R). 

47. A data processing system according to claim 46, wherein, in case the responses are 
absorbances, the absorbances being determined by calculating 



wherein the intensities (I) are obtained in a measuring region of the slave instrument by: 



an intensity (I), if necessary corrected with respect to dark current of the 





exposing an object in the riieasuring region to low and high X-ray energies 
and detecting with detectors the intensities Isampie(low) and l S am P !e(high) 
respectively 

detecting with the detectors the intensities l da rk(low) and l da rk(high) from said 
detectors when no radiation reaches them; 



and 
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exposing said detectors to the low and high X-ray energies signals when no 
object is present in the measuring region and detecting lair(low) and 
Uir(high). 

48. A correcting system comprising a slave instrument for obtaining responses and a 
data processing system for correcting the responses, the responses representing 
measurement performed with the slave instruments and the responses being based on 
detecting by the slave instrument attenuation and/or reflection and/or scatter of 
electromagnetic radiation, in particular X-ray, in/from a object exposed to said 
electromagnetic radiation in at least two spectral ranges, the set of responses comprises 
one or more pair of related responses (Q tow s and Qhigh*), said correcting system 
comprises 

processor means for determining the one or more pair of related responses 
(Qiow* and Q high s ) based on measurement on an object with the slave 
instrument, 

means comprising processor means adapted to perform a correction of 
responses by utilizing a correcting according to any of the claims 22-35, - 
said processor means comprises 

means for accessing memory means storing a first set of correction 
coefficients (B 0 ; B 2 ) 

processor means for determining the ratio [Q low / Qhigh] 00 "" by the correcting 
function; 

processor means for determining the corrected high energy response 
Qh>gh corr by the further correcting function; 

and 

processor means for determining the corrected low energy response CW 00 "" 
by multiplying [Q low / Q ww J 00,r by 0^°™. 

49. A system according to claim 48, wherein each of the responses (Q) is either 

an intensity (I), if necessary corrected with respect to dark current of the 
detectors; 

a transmittance (T) being derived from intensity as a ratio between intensity 
resulting from measuring on an object and reference intensity; 
an absorbance being defined as the negative logarithm to a transmittance 
(A— log(T)) such as logarithm base 10, the natural logarithm, or any other 
logarithmic function; 
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or 

a reflectance (R) expressing the reflectance from the surface of the object, 
the reflectance (R) is preferably linearized using the Kubelka-Munk 
transform (K/S=(1-R)/2R). 

50. A system according to any of the claims 44-49, further comprising storage means 
wherein a set of responses (CW" 1 and Qhigh") for a set of stable objects measured on a 
master instrument are stored and /or storage means wherein the first set of correction 
coefficients (B 0 ; B t ; B 2 ) and/or the further correcting function is stored. 

51 A dual X-ray instrument comprising a system according to any of the claims 44-50 
adapted to carry out a method according to any of the claims 1-35. 

52. A set comprising one or more stable objects for, or used during, carrying out a 
method according to any of the claims 1-21, each object being characterized by being 
composed by at least two different chemical compositions which are substantially stable 
and each stable object is having response, such as absorbance, properties which are - 
similar to the response, such as absorbance, properties of an object subjected to the 
method according to any of the claims 22-35. 



53. A set of stable objects according to claim 52, wherein for each of the stable object a 
first member of the at least two different chemical compositions is one having X-ray 
response properties, such as absorbance properties, similar to adipose tissue, and a 
second member of the at least two different chemical compositions is one having X-ray 
response, such as absorbance, properties similar to muscle tissue. 

54. A set of stable objects according to claims 52 or 53, comprising a plurality of stable 
objects have varying thickness and/or area! density. 

55. A set of stable objects according to claim 54, wherein the plurality of stable objects 
being integrated into a single stepped item. 



56. A set of stable objects according to any of the claims 52-55, wherein each object 
comprised in the set of objects is stable in the sense that the X-ray response, such as 
absorption, properties of the object does not change more than 0.1 %, such as no more 
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than 0.01 %, such as no more than.0,001%.. within at least 10 days, such as at least 1 
month, preferably at least 1 year. 

57. A plurality of stable objects according to any of the claims 52-56, wherein the number 
of stable objects in the set of stable objects are at least 8, such as at least 12, preferably 
at least 15, or even at least 20, such as at least 26. 
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Application of multivariate calibration models for the determination of the properties of the 
batch of meat: 

For all points: Calculate the fat content. 

Calculate the areal density. 

Multiply fat content and areal density, generating a "fat map". 
Add all points in the "fat map" to give the total fat weight, (Fbtoi) of the batch. 
Add all areal densities for the batch to give the total fat weight, (Wtotai ) of the batch. 
Calculate average fat content of batch as the ratio Ftotai/Wtotai 
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Application of multivariate a regression method (PLS) to the two set of 
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and application of linear regression to determine the coefficients a, 0 in Abigh corr 
= a * Ahigh s + 0. Store all coefficients in a memory in the computer in the X- 
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Application of multivariate calibration models for the determination of the properties of the 
batch of meat: \r-3J6 
For all points: Calculate the fat content. 

Calculate the areal density. 

Multiply fat content and areal density, generating a "fat map". 
Add all points in the "fat map" to give the total fat weight, (F to tai) of the batch. 
Add all areal densities for the batch to give the total fat weight, (Wtotai ) of the batch. 
Calculate average fat content of batch as the ratio Ftotai/Wtotat 
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